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A zero temperature dynamics of Ising spin glasses and ferromagnets on random graphs of finite 
connectivity is considered, like granular media these systems have an extensive entropy of metastable 
states. We consider the problem of what energy a randomly prepared spin system falls to before 
becoming stuck in a metastable state. We then introduce a tapping mechanism, analogous to that 
of real experiments on granular media, this tapping, corresponding to flipping simultaneously any 
spin with probability p, leads to stationary regime with a steady state energy E(p). We explicitly 
solve this problem for the one dimensional ferromagnet and ±J spin glass and carry out extensive 
numerical simulations for spin systems of higher connectivity. The link with the density of metastable 
states at fixed energy and the idea of Edwards that one may construct a thermodynamics with a 
flat measure over metastable states is discussed. In addition our simulations on the ferromagnetic 
systems reveal a novel first order transition, whereas the usual thermodynamic transition on these 
graphs is second order. 



PACS numbers: 05.20, 75.10 Nr, 81.05 Rm. 



I. INTRODUCTION 

Recently there has been much experimental and the- 
oretical interest in the properties of granular media. In 
such systems the thermal energy available is not suffi- 
cient to allow the rearrangement of a single particle and 
hence the system is effectively at zero temperature in the 
thermal sense. The fact that the problem is not trivial 
lies in the fact that such systems have an exponentially 
large number of such metastable states, which may be 
also called blocked or jammed configurations. Edwards 
associated an entropy to these configurations 



Ssdw = ln{N M s) 



(1) 



where Nms is the total number of metastable states of 
the system Q. It is reasonable to assume that in com- 
plex systems such as granular media SEdw is extensive 
meaning that Nms = exp(iVs) where s is the entropy 
per particle, alternatively one may work with an entropy 
per unit of volume which is clearly a more natural choice 
in granular media. Because the system has an extensive 
number of blocked configurations, if it is prepared from 
a random initial state it will lower its energy via only 
energy lowering rearrangements until it becomes stuck 
in a metastable state. Normally this first encountered 
blocked state will not be that of lowest energy (or most 
dense packing). In order to change the state of the sys- 
tem an external perturbation such as tapping or shearing 
is required. In between perturbations the system relaxes 
into new configurations. A natural and practically very 
important question concerning this sort of dynamics is: 
what are the properties of the steady state regime ob- 
tained via such mechanical perturbation schemes ? 

Recently it has been shown that spin glasses and ferro- 
magnets on random graphs have an extensive entropy of 



metastable states and one may calculate this entropy at 
fixed values of the energy ||||. Therefore, though they 
are quite different physically to granular materials, these 
systems have an extensive entropy of metastable states 
as do granular media. The motivation of this paper is to 
see if one can understand certain steady state properties 
of mechanically perturbed systems in terms of their or- 
ganization of metastable states. The possibility of using 
spin glasses as a paradigm for granular material was first 
introduced in 

Let us recall an example of an experiment on a sys- 
tem of hard spheres reported in ||. A system of dry 
hard soda glass spheres is placed in a glass tube. The 
system is tapped by using a piston to move the tube ver- 
tically through a sine cycle. The tapping parameter T is 
defined to be the ratio of the maximal acceleration due 
to the piston in the cycle to g the acceleration due to 
gravity. After an initial irreversible curve, obtained by 
increasing the tapping rate slowly, the system arrives on 
a reversible curve where the density is a monotonic func- 
tion of r, the highest packing densities being obtained at 
lowest tapping rate. Numerical simulations on granular 
media Q reveal similar behavior (though the irreversible 
part of the experimental curve corresponding to a loosely 
packed fluffy metastable state was not seen). It was also 
observed that at small tapping the relaxation to the final 
density is extremely slow and is well fitted by an inverse 
logarithmic decay of the form 



P(t) = Poo~ 



Ap 



l + Bba(l + t/r) 



(2) 



where p^ (the final density) , Ap^ , r (the characteristic 
relaxation time) and B are fitting parameters. It should 
however be remarked that the behavior of granular sys- 
tems is strongly dependent on the tapping mechanism 
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and that horizontal shearing M leads to behavior quali- 
tatively different to vertical tapping. 

In this paper we extend and elaborate a preliminary 
report of the results of ||]. The philosophy of the paper 
is to examine spin glasses as paradigms for granular me- 
dia. Here the quantity corresponding to the density is 
the energy of the system. We allow the system to evolve 
under a random sequential zero temperature single spin 
flip dynamics where only moves which reduce the energy 
are allowed. When the system is blocked we tap it with 
strength p £ [0, 1/2], that is to say each spin is nipped 
with a probability p, the updating at this point being 
parallel. The system is then evolved by the zero tem- 
perature dynamics until it becomes once again stuck, the 
tapping is then repeated. Physically this corresponds to 
assuming that in granular media the relaxation time to 
a new metastable state is much shorter than the time 
between taps. A similar, though not identical, tapping 
dynamics has also been introduced independently in the 
context of three spin ferromagnetic interactions on thin 
hypergraphs [|| , also in the goal of studying the dynam- 
ics of granular media. We find that a stationary regime is 
reached after a sufficiently large number of taps, charac- 
terized by a steady state energy E(p) (analogous to the 
stationary density - the same analogy as used in j^] ) . The 
initial dynamics from the random initial configuration 
into the first metastable state is examined analytically 
for the one dimensional ±J spin glass or ferromagnet 
(the two are equivalent by a gauge transformation). We 
call this the initial fall and the average energy of the first 
metastable state visited Ef is computed. We then de- 
velop a mean field theory for the dynamics under falling 
then tapping, interestingly this theory appears to be ex- 
act in the case of the one dimensional system and one 
may calculate E(p) within this scheme, the results being 
in excellent agreement with the numerical simulations. 

Numerically we examine the tapping of spin glasses 
and ferromagnets of higher connectivity. For the spin 
glass we find that E(p) is, as in the experiments, a de- 
creasing function of p. For small p we define the ex- 
ponent 9 by E(p) — E(0 + ) + Ap e ', with A constant. 
In the one dimensional case we show analytically that 
E{p) ~ — l + ^2p, hence 9 = 1/2, whereas for spin glasses 
on thin graphs for connectivity superior to two we find 
that 9 — 1. However for p < 0.05 we find that the time 
to reach the steady state is extremely long and not acces- 
sible numerically. In this slow dynamical regime we find 
a slow relaxation of the time dependent energy, reminis- 
cent of that observed in experiments on granular media 
|^| and hence compatible with Eq. (||). 

In the case of the ferromagnet we find numerically that 
there exists a critical value p c of p such that for p > p c , 
E(p) > Eqs where Eqs is the energy of the ground 
state and the inequality is strict, and that for p < p c 
E(p) = Eqs- Hence in the ferromagnetic system there 
is a first order phase transition under tapping dynamics 



(in contrast to the usual thermodynamic ferromagnetic 
transition in these systems which is second order ]To[|). 

There have of course been many models studied to 
understand the compaction process in granular media 
p"T| |l3| |, which reproduce many of the experimental fea- 
tures. Here the spin glass is clearly far from a realistic 
realization of a granular media, however the fact that 
it has extensive entropy of blocked states and the obvi- 
ously natural form of the tapping dynamics implemented 
makes it a natural testing ground for ideas about dy- 
namics and possible thermodynamics of systems such as 
granular media. Moreover, it has been argued in |Tl] ] 
that the slow compaction regime is well explained if we 
assume that particles can rearrange themselves in such 
a way as to create a particle size void, which is quickly 
filled by a new grain. This mechanism involves crossing 
of energy barriers and leads to a logarithmic compaction 
before the asymptotic steady state regime Q. We ex- 
pect that the local rearrangements which occur during 
the tapping dynamics on spin glasses random graphs will 
be lead to behavior analogous to the slow glassy dynam- 
ics of systems as granular media. 

Of course one would ultimately like to obtain a the- 
oretical understanding of the asymptotic, steady state 
regime of lightly tapped granular media. Edwards has 
proposed jjj that a light tapping dynamics on granular 
type systems leads to a steady state whose properties 
are determined by a flat measure over the blocked or 
metastable states satisfying the macroscopic constraints 
involved (e.g. fixed internal energy and compactivity) . 
This idea has recently attracted much interest and has 
been examined in the context of various models [|l4|-p"7[ . 
In this paper we shall concentrate simply on the asymp- 
totic energy of the final tapped state, the study of the 
dynamics leading to this final regime is deferred for fur- 
ther investigation Gaj. We shall see that the calculation 
of the Edwards entropy as a function of energy gives us 
a possible explanation of the first order ferromagnetic 
transition. 

II. SPIN SYSTEMS ON THIN GRAPHS 

The models we shall consider are spin systems on ran- 
dom thin graphs. A random thin graph is a collection 
of N points, each point being linked to exactly c of 
its neighbors, c therefore being the connectivity of the 
graph. The distribution of metastable states in these 
systems has been recently considered in |2|||. The spin 
glass/ferromagnet model we shall consider has the Hamil- 
tonian 

H = — — ./, , a, jS.S , (3) 

where the Si are Ising spins, n,j is equal to one if the 
sites i and j are connected and zero otherwise. The 
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fact that the local connectivity is fixed as c imposes 
the local constraints ^ ■ riy = c, for all sites i. In the 
spin glass case the Jy are taken from a binary distribu- 



tion where </;, 



-1 with probability half and = 1 



with probability half. In the ferromagnetic case Jij = 1. 
Here we define a metastable state as a spin configuration 
where any single spin flip does not increase the energy 
of the system. Mathematically the total number of these 
metastable states is expressed as: 



iV 



n ms = ti'[[o [ 



(4) 



It should be pointed out here that the definition of 
metastable states is of course dependent on the dynamics 
of the system, in contrast with micro-states in classical 
statistical mechanics. Whether or not, in certain cases, 
the information about the dynamics encoded in the cal- 
culation of the entropy of metastable states is enough to 
allow one to predict the properties of the steady state 
regime is an open question. 

The fact that, in our definition of metastable states, 
we include the marginal case (where the energy change is 
zero) implies that here 0(x), the Heaviside step function, 
is taken such that 8(0) = 1. In the context of granular 
media, where friction plays an important role, this is a 
natural choice as one certainly needs a non zero force in 
order to make a grain move. With this definition, the 
total number of metastable states of internal energy E 
per spin is formally given by 

N MS (E) = Tr J] 6 | JijrHjStSj j S(H - NE) (5) 
»=i J 

The corresponding Edwards entropy per spin, at fixed 
energy E per spin, is then given by 

s Edw (E)= In (N MS (E))/N (6) 



III. THE ONE DIMENSIONAL 
FERROMAGNET/SPIN GLASS 

We remark that by a gauge transformation the one di- 
mensional ferromagnet and ± J spin glass are equivalent 
and place ourselves for transparency in the context of 
the ferromagnet. Let us remark that the zero tempera- 
ture Glauber dynamics of the one dimensional ferromag- 
net can be explicitly solved Jl9| , here diffusion of domain 
walls occurs and the dynamics does not get blocked. In 
the Glauber case one may close the dynamical equations, 
however here such a closure scheme does not seem possi- 
ble. The zero temperature Kawasaki dynamics (conserv- 
ing the total magnetization) of the one dimensional Ising 



model, where the system can freeze, has been solved in 

To solve the dynamics of the one dimensional ferro- 
magnet we consider the dynamics from the point of view 
of the bonds. We define a fault of length n to be a se- 
quence of n neighboring adjacent domain walls. The zero 
temperature dynamics takes place within these faults via 
the flipping of one of the n — 1 spins contained between 
the n domain walls. We define Ii(ri) to be the indicator 
function that starting from bond i there are exactly n 
consecutive domain walls (there being no domain wall on 
bond i — 1 and no domain wall on bond n + i but all 
the intervening bonds have a domain wall). In the initial 
configuration we take the probability that a given spin is 
different to its left neighbor (that is to say the probabil- 
ity of a domain being between two spins) to be a. Hence 
if a = we have an initially ferromagnetic configuration, 
if a = 1 it is an antiferromagnetic configuration, the case 
a = 1/2 corresponds to a completely random configura- 
tion of maximal entropy. The total energy of the initial 
configuration is then given by 



So = -iY + 2]r^(n) 



(7) 



as the energy is given by the ground state energy — N 
plus two times the number of domain walls (excitations). 
Defining Kq(ti) to be the probability that one has n con- 
secutive domain walls starting from a given site (hence 
Ko(n) = (I(n)) where (•) indicates the average over the 
initial conditions), with the initial conditions introduced 
above one finds 



2„n 



K (n,N) = (1 -a) 2 a 



(8) 



(The length of the defaults has a geometric distribution) . 
The initial energy per site of a configuration generated in 
this manner is therefore (using the translational invari- 
ance of the system when N — > oo) given by 



En = -1 



9 N 

-Y 

n=l 



K (n, N)n 



(9) 



For the distribution of initial configurations considered 
here we find therefore that 



E 



-1 + 2a 



(10) 



We define by x( n ) the average number of isolated do- 
main walls left by a fault of n consecutive domain walls 
after the zero temperature dynamics described above has 
finished. The final energy of the system per site Ef is 
therefore 



2 N 

Ef = -l + ^J2 K ^ N )x(n) 

n=l 



(11) 
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By the definition of the spin dynamics, domain walls 
disappear by pairs of two neighboring domain walls. It 
is clear that x(l) = 1) x(2) = and we set x(0) = 0. 
Within such a fault the dynamics proceeds by flipping 
one of the n — 1 spins between the domain walls. By 
recurrence, after a random flip we obtain 



1 ™ _1 

X(n) = -J2x(k-l) + x(n-k-l) ( 

n — 1 z — ' 



121 



fc=i 



We solve equation (|T^) by introducing the generating 
functional 



The resulting equation for g(z) is 



dg 



2z 2 



1 



dz \(1 — z) 2 z 



+ - g(z) 



(13) 



(14) 



Solving this with the appropriate boundary conditions 
one obtains 



g[z) 



z exp(— 2z) 



and substituting this into Eq. (fyl) yields 



(15) 



E f = -1 + 2(1 - a) 2 X(n)a n = -1 + 2(1 - afg{ 



n=l 



thus giving the result 



Er 



-1 + 2aexp(— 2a) 



(16) 



(17) 



This yields a value of Et for the completely random 
initial configuration where, a — 1/2, of —0.632121. In 
fact the value of Ef is maximal for the case a = 1/2. 
For the totally antiferromagnetic initial condition, where 
= 1, here we find Ef — —0.72933. Clearly when a = 
the system is already in its ground state and we find 
Ef = — 1 as we should. We note that these values (and 
those for all a) have been checked with and are in perfect 
agreement with our numerical simulations. 

This calculation with the one dimensional fcrromagnet 
demonstrates two important points: 

• The final value of the energy Ef depends strongly 
on the initial configuration. 

• The system does not fall into a state of energy 
corresponding to the maximum of Nms(E). In 
HU it was shown that N MS (E) ~ exp(Ns(E)) 
where s(E) is a concave function peaked at E* = 
-1/V5 « 0.44721. Hence even if the total number 
of metastable states is dominated (in the statistical 



sense) by those of energy E* , generic initial condi- 
tions always seem to lead to an energy lower than 
this |Q. In [|l| the value of E f for a variety 
of zero temperature dynamics (sequential, greedy 
and reluctant) in the fully connected Sherrington 
Kirkpatrick (SK) spin glass model [^2| was stud- 
ied, similar behavior was found. 

When the one dimensional system is tapped we find re- 
sults in line with those described later for the spin glasses 
at higher connectivity. The curve of E(p), the asymptotic 
stationary value of the energy at a given p, is shown in 
Fig.(|l|) from 100 systems of size 10000 spins. 

The time taken to reach a stationary value for E(p) 
were rapid for larger p but for small p there is a very slow 
relaxation to the final asymptotic state which is of the 
form 1/y/i, where t is the number of taps. This is easily 
understood as at very slight tapping order p effects dom- 
inate at early times, this means that: (i) Isolated pairs 
of domain walls within large domains are immediately 
destroyed once tapping is stopped, (ii) Flipping a spin 
either side of a domain wall creates domain wall diffu- 
sion and with this annihilation by coalescence of domain 
walls. Hence the dynamics at small p and early times is 
qualitatively the same as that for the low temperature 
Ising model coarsening [23] ]. 

In order to go beyond our first calculation of Ef and 
solve the tapping dynamics we consider a mean field the- 
ory for the dynamics of a system of connectivity c. We 
shall see that at c = 2 this theory gives the analytic re- 
sult ( |l7|) and reproduces to within numerical errors the 
numerical tapping results. Once again we concentrate 
on the dynamics on the bonds. For a given site define 
x to be the difference between the number of unsatisfied 
and satisfied bonds. Hence x is the local field on the 
spin at this site and x G — c, — c + 2, • • • , c — 2, c. If 
x > then the spin can flip bringing about the change 
x — > —x. In addition we denote by P(x, k) the prob- 
ability that the site of interest has local field x after a 
total of k attempted random sequential spin flips under 
the zero temperature falling dynamics (that is to say the 
dynamics in between taps). We define /+ and /_ the 
probabilities that a given spin can flip conditional on the 
fact that the bond with a given neighboring site is not 
satisfied or satisfied respectively. Formally we have 

f± = Prob(x > 1 given bond is not satisfied/satisfied) 

(18) 

We may turn around this conditional probability using 
Bayes' Theorem to obtain 



f± = 



Prob (x > and given bond is not satisfied/satisfied) 
Prob( given bond is not satisfied/satisfied) 

(19) 
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Given that a site has local field x, it must have (c + x)/2 
unsatisfied bonds and (c—x)/2 satisfied bonds. Therefore 
we find that 

Prob (x > and given bond is not satisfied/satisfied) 
= ^P(x)(c±x)/2c 



x>0 



and 

Prob (given bond is not satisfied/satisfied) 
= ^2P{x)(c±x)/2c 

X 

Putting these results into Eq. ( |l8| ) then gives 
E* >0 p (*)(c±z) 



f± = 



E a ^)(c±. 



(20) 



If we are interested in the spin at site i the possibilities 
between time k and k + 1 are 

• The spin at site i is chosen and x > 0, then the 
spin at site i will flip and x goes to —x. 

• The spin at site i is chosen and x < 0, then the 
spin at site i can not flip and x does not change. 

• A neighbor of site i with positive local field is cho- 
sen and so flips. In this case, x goes to x + 2 or to 
x — 2 depending whether or not the bond with site 
i was satisfied or not satisfied. 

• A neighbor of site i with negative or zero local field 
is chosen and so does not flip. In this Ccise, x does 
not change 

• One chooses neither the spin at site i nor any of its 
neighbors, and so x stays x. 

Assuming that the distribution at every site is given by 
P(x, k) and assuming independence between the values 
of x from site to site (the mean field approximation) we 
obtain 

P(x, k+l)= °-^P(-x, k) + d -^-P(x, k) 



N 

N-c-1 
N 



N 



P(x,k) 



^V.A-)(^(l-/ +) + ^(l-/_) 



P( x + 2,k) C -±^f + 
P(*-2,k) C -^f_ 



(21) 



In this equation we have to define 8(0), the choice com- 
patible with the conservation of probability is 8(0) = 1/2. 
Taking the limit N — > oo we may introduce the continu- 
ous time r = k/N and obtain 



= 8(~x)P(-x) + 8(-x)P(x) - (c + l)P(x) (22) 

ar 

+ p(x) (£±£(i-/ + ) + £z£(i_/_) 

+ P{x + 2) C -±^f + + P(x - 2) £^±2/_ 



The average energy per site at time r is then given by 
E(t) = i Ex xP(x, t) and one can show that the above 
mean field equation (^||) respects the exact identity for 
the evolution of the average energy per spin 



dE 



-2^>P(x,t). 



(23) 



x>0 



The case where c = 2 (the one dimensional case) is ac- 
cessible to analytic solution and we proceed by defining 

u(t) = P(-2,t) 
v(t) = P(0,t) 
w(t) = P(2,t) 

One finds that 



f-=0 



(24) 



2w 



v + 2w 

and the full mean field evolution equations become 
da 



dr 
dv 

dw 
d^ 



f+v + w 
-f+v + 2f+w 
-(2f + + l)w 



(25) 



If we look for a stationary solution of (|2J) and (|25|), 
we find w = 0, which expresses the fact that when r is 
infinite, the system is in a metastable state. To solve (|2~i| ) 
and (|25|), we introduce A so that v — Xw. Then A obeys 
the equation: 



dX 
d7 



A + 2, 



(26) 



with the initial condition A(0) = v(0)/w(0). Then (g|) 
becomes: 



du 2Xw 



w 



dr A + 2 

dw A + 6 

— = w 

dr A + 2 

A + 2 = (A(0) + 2) e T 
and w(t) and v(t) are given by 



(27) 



(28) 
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w(t) = w(0) exp U + ^^(^ r - 1)) (29) 
v{t) — —2w(r) 

+ w(0)(X(0) + 2) exp ( A(Q) 4 +2 (e~ T - l)j (30) 

The probability to have a positive value for the local fields 
then goes to zero at infinite r as expected and the limit 

4 

of v is u(oo) = (u(0) + 2w(0)) e A <°)+ 2 . If we consider 
the geometric initial conditions used in the previous ex- 
act calculation of Ef, the induced initial conditions are: 
u(0) = (1 - a) 2 , v{0) = 2a(l - a) and to(0) = a 2 . In 
this case we obtain Ef = — I + u(oo) = — I + 2ae~ 2a 
reproducing the exact result (^)- Tapping the sys- 
tem with tapping probability p, starting from the values 
{u(oo), v(oo), w(oo)}, we obtain the new tapped values 
{t/(0),i/(0),«/(0)}. Defining q = (1 - p), the relations 
between the old and tapped probabilities are: 

u'(0) = (l-3pq)u(oo)+pqv(oo) (31) 
v'(0) = 2pq u(oo) + (1 - 2p<?) u(oo) 
tf/(0) = 

Then, after another zero temperature evolution of the 
system, it reaches a new local energy probability distri- 
bution with w'(oo) = and: 

u'(oo) = (Apq + v(oo)(l — ipq)) 

x -P ('I r% 3 ^ (32) 

F V 4pq + i;(oo)(l-4pq)y V 7 

At this stage of computation one should remark that 
this recursive equation contains one of the main features 
of our numerical simulations, that is reversibility. In- 
deed, the process involved in (|3^) will reach an asymp- 
totic value which is independent of the initial conditions 
and depends on p. Hence in the steady state regime under 
tapping, the probability v s (p) (the subscript s indicating 
steady state) for sites to have zero local field is solution 
of the fixed-point equation: 

Vs(p) = (4pq + v s (p)(l- Apq)) 

x exp (— p^, —) (33) 

F V ipq + V s (p)(l-4pq)J 

This equation can be solved numerically and the result 
is shown in Fig. (p]) in comparison with the numerical 
simulations which we see is excellent. 

The small p behavior of E(p) from ( ^3| ) is: E(p) = 
— 1 + ^/2p+ 0(p), indicating that in this case 9 = 1/2. 

Given the mean field nature of the above calculation we 
have used we do not expect this approximation to cor- 
rectly describe the approach towards the steady state, 
by direct comparison with the numerical simulations we 
have verified that this is indeed the case. Let us remark 
here that a defect of the mean field approximation scheme 



is that it cannot distinguish between a spin glass and a 
ferromagnet, this is clearly not a problem for the one 
dimensional situation where the two are identical. 



IV. HIGHER CONNECTIVITIES 

The systems which we study are ± J spin glasses or uni- 
form ferromagnets on random graphs with fixed connec- 
tivity c. Let us first recall some analytical results of fl[| . 
It has been found that the mean number of metastable 
states increases exponentially with the number of sites in 
both cases. In addition in || an annealed approximation 
to the Edwards entropy per spin of metastable states at 
fixed energy E was carried out: 

s E dw{E) = ]n({N MS {E))) /N (34) 

which may be exact for the ferromagnet as the calculated 
entropy is always positive. Moreover, there is an energy 
threshold E* above which the results are the same for the 
±J spin glass and the ferromagnet; below the ferromag- 
net has more metastable states and a non-zero magneti- 
zation. Hence, as far as the energy density of metastable 
states is concerned, both ferromagnet and spin glass are 
the same above E* - that is the effect of loop frustration 
is negligible. In this regime, one also suspects that the 
zero temperature dynamics are the same. In particular, 
numerical simulations with 100 samples of TV = 10000 
sites for connectivities of 3, 4 and 5 have found the same 
Ef for the spin glass and ferromagnet with very good 
accuracy (the relative error is about 10~ 6 ). The results 
are show in Table (1). 

The result of tapping experiments on the systems with 
c = 3 is displayed in Fig.(g). There is some critical tap- 
ping rate, p c , above which the curves of E(p) versus p 
are the same for the spin glass and the ferromagnet. 
Moreover, the ferromagnet is subject to a phase tran- 
sition under tapping dynamics at p c such that for p < p c , 
the steady state reached is the ground state. Finite size 
effects have been studied and revealed that the transi- 
tion is first order (in as far that E FM (p+) ^ E FM (p-)), 
in contrast to the usual thermodynamic ferromagnetic 
transition in these systems which is second order [fl0|| . 
For the ferromagnet in the region close to p c one finds 
an excellent scaling of the energy as a function of N, 
E FM (N,p) — f(N(p — pc)) as shown in the inset of Fig. 

. This scaling may be used to optimize the determina- 
tion of p c . The first order nature of the transition may be 
seen explicitly by looking at the histogram over time (in 
the steady state regime) of the average energy per spin 
at p = p c ; one sees in Fig. (^) two separated peaks in 
the distribution and not a single peak which splits into 
two as one would expect for a second order transition. 
Near p Cl for systems of finite size, there is therefore co- 
existence of the two phases. The time dependence of the 
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average energy per spin in the simulation leading to the 
histogram Fig. (^) is shown in Fig.(||). One sees that the 
system tunnels between the two coexisting states. The 
typical time for this tunneling increases as the system 
size increases, indicating, in thermodynamic language, a 
free energy barrier between the two phases. As the sys- 
tem size is increased the occupation of the intermediate 
states of energy between the two phase of energy Eqs an d 
E(pt) (between the two peaks in Fig.(||)) is suppressed. 

We have also measured E{p) by studying single sys- 
tems of very large size {N = 10 6 ), the results are shown 
in Fig. (H). Here again above p c the curve for the spin 
glass and the ferromagnet are completely indistinguish- 
able and the ferromagnet reaches the ground state below 
p c . For such large sizes, one no longer sees a coexistence 
of two phases around p c as presumably the tunneling time 
has become much larger than the simulation time. 

Moreover, for N — 10 6 , for p > p c the full temporal 
plots (and not just the steady state values) of E SG {p,t) 
and E FM (p, t) (where t is the number of taps) are indis- 
tinguishable. For p < p c the two curves are identical up 
till a time tdif, which depends on the initial configura- 
tion and the sequence of spins flipped during the tapping 
process, and diverge after tdif, when the ferromagnetic 
system reaches quickly the ground state (see Fig.(||)). 
Once the ferromagnetic system has broken the Zi sym- 
metry the easiest way to lower the energy is to flip the 
spins which are opposed to the global magnetization (be- 
cause they are more probable not to be in the direction 
of their local field) until all the spins are —1 or +1. 

Identical behavior was found in the systems with c = 4 
and 5. The comparison of E SG (p) and E FM (p) for c = 4 
is shown in Fig. (^|). Remark that if we compare the 
different values of p c when increasing c, and considering 
only odd (or even) connectivities, we find that p c grows, 
and we expect that it goes to 1/2 when c is very large, 
as the metastable states are more and more magnetized 
when c grows (in the case of the fully connected ferro- 
magnetic Ising model, only the two ground states are 
metastable, so p c — 1/2). 

The behavior of the spin glass systems is similar to 
that for the system with c = 1. The steady state en- 
ergy E SG (p) is a monotonically decreasing and contin- 
uous function of p. For small p one finds that here 
E SG (p) ~ E SG (0) + Ap giving 6 = 1 in contrast with 
9 = 1/2 in the one dimensional case. 

A tentative explanation for the ferromagnetic transi- 
tion: In H it was also shown that for the ferromagnet 
the Edwards entropy as a function of E is concave for 
E > E* and convex for E < E* . The value of E(p+) 
obtained from the tapping experiments are very close to 
those obtained for E* in ||], the energy at which s^dw 
becomes convex. The results are shown in Table (1). En- 
couraged by this striking observation we will try to make 
a tentative link with a possible thermodynamics for such 
systems. If we imagine that the energy of the system 



is governed by a partition function inspired by the flat 
Edwards measure over metastable states fp]|lF 



Z = J dEN MS (E) cxp(-N(3E) 



(35) 



where (3 is a Lagrange multiplier corresponding to the 
inverse Edwards temperature which depends solely on p 
and not on E and is a monotonically decreasing function 
of p for p 6 [0 + ,l/2]. The monotonicity hypothesis is 
supported by the simulation results that E(p) decreases 
with decreasing p. Clearly the energy which dominates 
in the sum is that obeying 



dsEdw(E) 
dE 



(3 = 



(36) 



However if this saddle gives a true maximum of the action 
one must have also that 



d 2 s Edw {E) 
dE 2 



< 



(37) 



and hence the Edwards entropy must be concave for 
the energy considered to be thermodynamically stable. 
Hence for E < E* , this suggests that the only stable 
energy is the ground state. 

If one would like to push the analogy of the thermo- 
dynamics of first order phase transitions to its limits and 
assume that (3 is a continuous function of p (it is clear 
that if this is not the case one could trivially obtain a 
first order transition), one would also expect that the 
free energy of the two phases is equal at f3 c = (3(p c ) and 
hence 

s E dw(E*) - [3(p c )E* = s Edw {E GS ) - (3{ Pc )E GS (38) 

where Eqs is the energy of the ground state of the ferro- 
magnet given here by E G s = — c/2. It is also clear that 
SEdw(EGs) = hence one obtains 



P(Pc) 



SEdw(E*) 



E* — Eqs 

However one also has that 

PiPcj = P(E ) = — \e=e* 



and hence one should obtain 



dsEdw(E ) , 
dE 



\E=E* 



dE 



SEdw{E*) 



E* - E, 



GS 



(39) 



(40) 



(41) 



However the fact that the calculated annealed approxi- 
mation for SEdw(E) is convex for E < E* means that 



dsEdw(E) SEdw(E*) 
\E=E' > 



dE 



E* - E, 



GS 



(42) 



and hence the equality (^) is not respected. It is possi- 
ble that the exact quenched calculation of SEdw{E) would 
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give a different value from the annealed calculation of |i| , 
however one should not expect the result to be too dif- 
ferent qualitatively to that of 0. 

In addition we remark that the annealed calculation 
In ((Nms{E)}) gives an upper bound for the quenched 
value (In (N ms (E))} (from Jensen's inequality) . Also the 
value of E* found with this approximation seems to be 
greater than that found numerically for E(p^) by a very 
small amount (« .005 for c = 3). If we accept the pos- 
sibility that E(p+) = E*, a more probable hypothesis 
is the collapse of validity of Edwards' hypothesis in the 
region of energy where the Z2 symmetry is broken on 
metastable states, as has been found for the three di- 
mensional random field Ising model in [fl4|| . 

Let us mention here that we measured the energy in 
the simulations over a few hundred time steps after the 
energy appeared to stop to decay. To be sure that the 
systems considered here were in a stationary regime (and 
that the energy was not decaying extremely slowly as a 
function of time) we measured the correlation function 
at different waiting times t w (the number of taps after 
the initial preparation of the system) , that is to say 

1 - 

C(t + t w ,t w ) = —J2 s i( t n,)S i (t w +t) (43) 
i=i 

In the stationary regime this should be a function only of 
t. In out of equilibrium systems the fact that the system 
is not in equilibrium shows up strongly as aging in the 
correlation function i.e. C(t+t w ,t w ) ^ C(t) (see p4j and 
references within), even though the energy may be decay- 
ing so slowly that it appears to have reached its asymp- 
totic equilibrium value. The time translational invariance 
of C(t + t w , t w ) is thus quite a rigorous test of whether 
the steady state regime has been attained. For example 
for the case p = 0.02, with the waiting times t w — 30000 
and t w = 60000 is shown in Fig. (0), one sees, clearly 
that after the appropriate translation of the t axis, the 
two functions collapse perfectly onto one another. One 
also sees that the decay of C(t) (in the longtime regime 
we can now eliminate the t w dependence) is exponential 
at large t and also that C(t) decays to zero, indicating a 
form of ergodicity in the system. As pointed out in [pr[ , 
this behavior of the correlation function seems a neces- 
sary condition for the validity of the scenario of Edwards, 
that under tapping all metastable states satisfying the 
relevant macroscopic constraints (fixed energy and com- 
pactivity) are equiprobable in the stationary regime of 
gently tapped or perturbed system. 

In order to test the accuracy of the mean field approx- 
imation at high energies (where the system does not dis- 
tinguish between the ferromagnet and the spin glass), we 
have compared the value of Ef obtained in the numerical 
simulations with the result £7jr obtained by numerical 
integration of Eq. ( ^2|). The comparison is in Table (1) 
and we see that the agreement is quite good. 



Finally we mention that we have also examined the 
reversibility of the tapping mechanism. If the system 
is tapped for a sufficiently long time, compatible with 
the relaxation times discussed above, the system is com- 
pletely reversible. This reversibility was found in the ex- 
periments in H once the system had left the initial fluffy 
state. 

V. CONCLUSION 

Granular media are a natural example of systems hav- 
ing an extensive entropy of metastable states. In such 
systems the role of thermal fluctuations are negligible and 
in order to evolve one must apply some external tapping 
mechanism. One would ultimately like to be able to for- 
mulate some sort of thermodynamics for such systems. 
The proposition of Edwards (l) for a thermodynamics of 
such systems is an important step in this direction and 
has had some succe ss [|14|Jl7|| but it has been shown not to 
be generically true ]14|. A more general understanding of 
the asymptotic states of tapped systems has far reaching 
implications for computer science as the tapping mecha- 
nism studied here is similar to certain algorithms used in 
optimization problems. 

We have presented what appears to be an exact cal- 
culation of the steady state energy of a tapped one di- 
mensional spin glass or ferromagnet. For this problem 
we have obtained the fixed point equations for the dis- 
tribution of local fields under tapping. These equations 
also explain the reversibility observed in the numerical 
simulations. In a wide context of models we confirm the 
observations of |^,^,^|, that if one reduces the strength 
of tapping, then the compaction process, corresponding 
here to the reduction of the energy of the system, be- 
comes more efficient. The existence of a first order type 
phase transition for tapped ferromagnets on random thin 
graphs is of great interest, the possible explanation using 
the calculations of [|| on the Edwards entropy for this 
system indicates the possibility that one may eventually 
construct a more general theory for the thermodynam- 
ics and even phase transitions in tapped systems. One is 
tempted to speculate that generically the convexity of the 
Edwards entropy below a certain energy threshold (de- 
noted here by E*) leads to a collapse to the ground state 
energy Eas, the metastable states in the intervening en- 
ergy values being unable to support a stable thermody- 
namics. In terms of granular media this sort of transition 
would correspond to a transition between a random close 
packed state to a crystalline close packed state. It would 
be interesting to find other systems (both theoretical and 
experimental) showing the same collapse phenomena in 
order to test this idea. 

Finally let us mention some of the open questions 
posed by this study we believe to be of interest for fu- 
ture investigation. Clearly a general goal would be, in 



8 



the spirit of Edwards, to develop a thermodynamics to 
describe the stationary regime of tapped systems such as 
those studied here. The exact results presented on one 
dimensional systems here provide a completely analytic 
understanding of the tapping dynamics which one may be 
able to rederive from static considerations. Indeed it has 
been shown JlTj in this simple context that several steady 
state observables may be predicted using Edwards' mea- 
sure. The phase transition found in the case of the fer- 
romagnetic systems studied here is extremely novel, an 
analytic understanding of this phenomena would be de- 
sirable, perhaps there exists a percolation type argument 
which would allow one to evaluate p c . Also of interest 
is the decay of a system towards its final steady state 
energy. The slow logarithmic decay described by Eq. (||) 
has been used successfully to fit the experimental data 
of H and the simulation data of [PJpj. Our prelim- 
inary study of finite connectivity spin glasses and the 
SK model [H indicates the presence of a slow dynam- 
ical regime for small values of the tapping parameter p 
which is also compatible with a slow logarithmic decay, 
but the curves can also be well fitted by power law decays 
(with the same number of fitting parameters). There ex- 
ist however phenomenological arguments [[Lll and exact 
calculations and simulations on toy models 
support logarithmic decay. 
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Figure Captions 

Tab. 1. Comparison of the numerical values of E* and E(p+) for different values of the local connectivity c. The 
result for E* when c = 3 is a truncation of the analytical value —15/14. 

Fig. 1. Comparison between numerical simulations of tapping experiments (b) and the analytical result (a) obtained 
with ©. 

Fig. 2. Numerical simulations of tapping experiments for the spin glass (c) and the ferromagnet ((a), (b) and 
(d)) for c = 3 for N = 1000 ((c) and (d)), N = 2000 (b) and N = 10000 (a). The inset shows the scaling 
E FM (N,p) = f(N(p - p c )) for p ~ p c for N = 400, N = 1000 and N = 2000. 

Fig. 3. Histogram of the energy per spin obtained from Fig. (Q) . 

Fig. 4. Single run for a ferromagnet of local connectivity c=3, for N = 10000 at p = 0.249 w p c . One sees that 
because the size is not too large, the energy switches between two values, one not far from that of the ground state 
Eqs and the other not far from E* . 

Fig. 5. Numerical simulations of tapping experiments for c = 3 ((c): spin glass, (d): ferromagnet) and c = 4 ((a): 
spin glass, (b): ferromagnet). Here, we have computed the asymptotic energy for only one sample of very big size 
N = 10 6 . The results are quite the same as for N = 10 5 and do not change if we average over several samples, which 
indicates that at these sizes, we are very near the thermodynamic limit and we can study only single runs to compute 
the energy. 

Fig. 6. Comparison of the energy versus time (number of taps) for the ± J spin glass (a) and the ferromagnet (b) 
for N — 10 6 spins at p = 0.05. We have displayed the magnetization for the ferromagnet (c), whose absolute value 



9 



increases with time, whereas that of the spin glass remains zero. 

Fig. 7. Correlation function for the ± J spin-glass with local connectivity c = 3 versus number of taps t for two 
values of the waiting time t w : t w — 30000 and t w = 60000. The tapping value is p — .02, the system contains 
TV = 1000 spins and we have averaged over 1000 samples. With triangles is shown the right part of the curve, which 
corresponds to t w = 60000, shifted to the left by 30000, it superimposes perfectly over the curve for t w = 30000 
demonstrating the time translation invariance of the correlation function. 
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